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Abstract 

A new nonlinear equation governing asymptotic dynamics of ripples 
is derived by using a short wave perturbative expansion on a generalized 
version of the Green-Naghdi system. It admits peakon solutions with 
amplitude, velocity and width in interrelation and static compacton so- 
lutions with amplitude and width in interrelation. Short wave pattern 
formation is shown to result from a balance between linear dispersion and 
nonlinearity. 



1 Introduction 

Since the classical works of Boussinesq [Q and Korteweg and de VriesQ, nonlin- 
ear evolution of long waves of small amplitude in shallow fluids has been widely 
studied. The asymptotic dynamics (in space and time) is by now well under- 
stood and represented by a large number of model equations, among which the 
different versions of the Boussinesq equations || , the Korteweg-de Vries (KdV) , 
the Benjamin-Bona-Mahoney-Peregrini (BBMP) Q, and the Camassa-Holm Q 
equations. 

On the contrary, short waves have been studied very little and only a few 
results are known on their asymptotic dynamics. The main purpose of this paper 
is to study nonlinear short surface waves in fluids, the ripples, shown to build 
up as a result of superposition of two surface motions: an oscillatory flow and 
a laminar flow. The oscillatory flow corresponds to mechanical perturbations 
which propagate like a wave. The laminar flow may be created in many ways: 
by the action of an external wind, or in a two-layer liquid where the upper one 
displaces particles belonging to the upper surface of the lower fluid, or by an 
external electric field if the surface particles are charged, etc... 
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Here we study an ideal fluid (inviscid, incompressible and without surface 
tension) in which surface displacement can be achieved through the action of a 
steady external wind directed parallel to the water surface. A surface wind on 
a lake which produces surface flow is an ideal physical environment conducive 
to the above mentioned phenomenon. 



2 Model equation for ripples propagation 

Let's consider u(x, t) which represents at time t an unidirectional surface wave 
propagation in the x direction of a fluid medium which is involved in a flow on 
a large scale. This large scale flow is a superposed motion of the surface which 
moves under the action of the wind with a velocity Cq in relation to the bulk. 
Using perturbativc methods we will show that ripples can propagate and obey 
the nonlinear equation 

u xt = --p-u- uu xx + {u x ) 2 (2.1) 
hc 

where subscripts denote partial derivatives, h is the unperturbed initial depth 



and g is the acceleration of gravity. The equation (2.1) has several types of 
interesting solutions. The peakon solution is (using the formula (1 — d zz )e~\ z \ = 
26(z), |) 

u = -aA 2 exp(-| ), a = (2.2) 

A hco 

where the width A is a free parameter. Unlike the peakon of Camassa-Holm 
equation, the amplitude (— aA 2 ), the velocity (aA 2 ) and the width (A) are in- 
terrelated. 

The static compacton solution is 

u= -8aA 2 cos 2 (4), |vI<2tt, (2.3) 
4A A 

and u = otherwise. Unlike the compacton solution recently introduced and 
investigated by Roseneau and Hyman B, this solution presents a dependence 



between width and amplitude. Solutions (2.2) is a coherent struture analogous 
to the solitary wave of KdV. Moreover, contrary to KdV, (2.1) possess a plane 
monochromatic wave solution of arbitrary amplitude A 

u = A expi(kx - Qt) (2.4) 

with the dispersion relation fl(k) — 3g/khco identical to the one of the linearized 
system . Moreover this dispersion relation is a function of k which has a good 
behavior in the short wave limit. 



Despite the fact that solution ( |2.3| ) is unnatural in the physics under consid- 
eration, it shows that ( |2 . l| ) is an adequate mathematical tool to modelize statics 
patterns in nature. 
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The underlying mechanism responsible for structural stability of solutions 
of the Camassa-Holm equation or the KdV equation with nonlinear dispersion 
is the balance between nonlinear dispersion, nonlinear convection and nonlin- 
earity. Indeed these equations are nonlinear evolution equations without linear 
dispersion (the plane wave is not a solution of the linear associated evolution 
equations). However the exibited solutions of (2.1) come from the balance be- 
tween linear dispersion and t wo n onlinear terms (in the case of ( |2.2[ )) and only 
nonlincarity (in the case of ( |2.3|) ). An essential point, which remains to be 
proved by numerical or analytical methods is: under what conditions do the 
solution ( pyi ) dominates the initial value problem of equa tion (2.1)?. Another 
important open problem is the structural stability of ( |2.2| ) and (2.3) . Here we 
only show that the plane wave ( |2.4|) is unstable. 



3 Physical context: modified Green-Naghdi sys- 
tem of equations. 

We consider an inviscid, incompressible, homogeneous fluid with density a . Let 
the particles of this continuum medium be identified by a fixed rectangular 
Cartesian system of center O and axes (x±, £2, £3) = (x,z,y) with Oz the up- 
ward vertical direction. We assume symmetry in y and we will only consider 
a sheet of fluid in the xz plane. This fluid sheet is moving in a domain with 
a rigid botton at z — and an upper free surface at z = <p(x,t). The vector 
velocity is v = (vi,«2) = (u,w) for < z < <f>{x,t). Thanks to homogenity and 
incompressibility the continuity equation reduces to 

V • v = u x (x,z,t) + w z (x, z,t) = 0, < z < 4>(x,t), (3-1) 

where V = (d x ,d z ). The Euler equations of motions (law of conservation of 
momentun) of a fluid under gravity g and for < z < (f>(x,t) are 

au(x, z, t) = — p x (x,z,t), (3-2) 

aw(x, z,t) — —p* z (x, z,t) — ga, (3-3) 

where p*(x,z,t) is the pressure and a superposed dot denotes the material 
derivative for x, z fixed. 

We complete now the fundamental equations of continuity and momentun 
conservation with appropiate kinematic and dynamics boundaries conditions. 
Let S(x, z, t) be the interface between the inviscid fluid sheet and the air (ex- 
ternal medium). We represented S(x, z,t) by the (classical) equation 

S(x,z,t) = (j>{x,t) - z = 0. (3.4) 

The kinematic condition is that the normal velocity of the surface S(x, z, t) must 
be equal to the velocity of the fluid sheet normal to the surface. The normal 
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velocity of the surface is 



S t 



while the velocity on the surface z = <p(x, t) is 

v = (u + Co, w), z = 4>(x, t) , 
whose normal component is 

II II 



From (3.5), (3.6) and ( p . 7| ) the kinematic conditions read 

4> t + u<t> x - w + c 4> x = o, z = 4>. 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



The equation ( p.6| ) leading to ( |3.8| ) lies at the heart of our approach. It points 
out that an external agent - wind in our case - drives the particles belonging to 
S(x,z,t). The motion is uniform, of value cq and in the x direction only. So 
S(x, Z, t) is a surface of discontinuity for the u component of v, which experiences 
finite jumps of value Co in z = <f>. This wind induced motion of the fluid's surface 
is found in particular in the hydrodynamics of lakes || . From a theoretical point 
of view this phenomenon was studied in the ray tracing theory by Lighthill in 
[ fio| where an equation analogous to (3.8) was derived. 

At the surface of the sheet z — (f>, there is a constant normal pressure po- 
At the bed z = 0, there is an unknown pressure p*(x,0,t) and the normal 
fluid velocity is zero: w = 0. Several types of approximations can be used 
in order to solve this wave problem. Here we adopt an approximation in the 
velocity field introduced by Green, Laws and Naghdi in and based on the 
monumental work of Naghdi |l^| . We assume that u is independent of z. This 
is equivalent to considering the vertival component w as a linear function of z. 
This simple and realistic assumption enables us to satisfy exactly the equation 
of incompressibility and the boundary condition at the bed. Hence u = u{x, t) 
and from (3.1) we have 



= z£(x, t) , £(x, t) = -u x (x, t). 



(3.9) 



Now we integrate (3.2) in the variable z from z — to z — <j) (the integration is 
granted by the Riemann's condition of integrability) . The result is 



a(u t + uu x )4> = -p x , 



where we use u = u t + uu x and 



p(x,t)= p*(x,z,t)dz - p (j>(x,t). 

Jo 



(3.10) 



(3.11) 
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Next, we multiply equation (3.3) by z and integrate from z = to z — (f> which 
yields 



(3.12) 



The pressure p can be eliminated using (3.12) in (3.10) and, with the help of 
y&M) we eventually obtain 



Ut + uu x + i 



(t>4>x{u xt + UU XX - (u x ) 2 ) + 
03 

-^-( u 3;t + l" 1 !! — ( u z) )a: 



(3.13) 



The remaining upper boundary condition on <j> reads using (3.9) 

<t>t + {<pu) x + c o 0x = 0. (3.14) 



With Co = 0, the system ( 3.13 ) ( 3.14 ) is the Green-Naghdi system of equa- 
tions fnl. The inclusion of the term cq(/) x drastically changes its dynamics: c$ 
cannot be eliminated neither by a Galilean transformation nor by a rescaling of 
u or <j). These extended Green-Naghdi equations represent the nonlinear inter- 
action between two separate forms of motion: a wave motion associated with 
the elastic response of the fluid to a perturbation, a surface (uniform) motion 
generated by an external agent. 



4 Nonlinear dynamics of ripples. 



Let u s consider now asymptotic nonlinear dynamics of ripples in ( 3.13|) and 
■ Ripples are, from a geometrical point of view, nearly local objects and we 
are looking for their large time behavior. So we need to introduce two appropiate 
variables: one space variable £ describing a local and asymptotic pattern and 
a time variable r measuring asymptotic time dynamics. These asymptotic are 
worked out by means of the change of variables 



C = -x 



et. 



(4.1) 



where the small parameter e is related to the size of the wavelength: £ — 2ir/k 



Such variables cannot always be defined as shown in 15 , but exist if the 



linear dispersion relation f2(fc) can be expanded as a Laurent series with a simple 
pole for k — > oo, that is 



n(k) = Ak 



D 

T 



c_ 



where A, B, C, are constants. Consequently the phase and group velocities 
remain bounded in the short wave limit k — > oo. This is the case here because the 
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dispersion relation for ( 3.1 3| ) (3.14), obtained for (j>{x, t) — h+aeiq>i[kx — Q(k)t] 
and u(x,t) — bexpi[kx — Q(k)t], has the asymptotic expansion 



= -<|k>i 



Co 



h 2 y h 



,3g 3g 2 . 1 1 . 

+ ^p+^- 



(4.2) 



Such an expansion justifies the change of variables (4.1) and it is essential in 
revealing in ( 3.13 ) ( 3.14 ) the asymptotic short wave dynamics fusf . We are 
looking now for nonlinear dynamics of ripples of small amplitude represented 
by the expansions 



= h + e 2 (H a + e 2 H 2 + e A H A + ...) , 
u = e 2 (U + e 2 U 2 + e i U 4 + ...) , 

which from ( 3.13 ) and ( |3.14 ) give at lower order in e: 

h(U ) C + Co{H ) c = , 

s(#o)c = y {(^o)rc + Uo(U ) cc - (C7b)|k 



(4.3) 
(4.4) 



(4.5) 
(4.6) 



Using (4.5) in ( |4.6| ), going back to the laboratory fields (</>, u) and co-ordinates 
(x, t) and integrating once we eventually arrive at equation (p.l|). 



5 Benjamin-Feir instability of the plane wave. 



For a plane wave, the nonlinears terms in (2.1) give rise to harmonics of the 
fundamental. Assume that a disturbance is present consisting of modes with 
sideband frequences and wavenumbers close to the fundamental. We can have 
interaction between harmonics and these sideband modes. This interaction is 
likely to produce a resonant phenomenon manifesting itself by the modulation 
of the plane wave solution. The exponencial growth in time of the modula- 
tion, originating from synchronous resonance between harmonics and sideband 
modes, leads to the Benjamin-Feir instability Jl^]. A formal solutions can be 
given via an asymptotic expansion conducing to the nonlinear Schodinger equa- 
tion (NLS) [O] . The particular interest of NLS is the existence of a general and 
simple criterion enable to detect stability or unstability of the monochromatic 



wave train. Let us seek for a solution of (2.1) under the form of a Fourier ex- 
pansion in harmonics of the fundamental expi(kx — fit) and where the Fourier 
components are developed in a Taylor serie in powers of a small parameter 7 
mesuring the amplitude of the fundamental 

i=p 00 
i=—p p=i 



In (5.1), = Ui ("star" denotes complex conjugation) and £ and r are 



slow variables introduced through the stretching £ = ^{x — vt) and r = 7 t 
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and where v will be determined as a solvability condition. The expansions 
( |5~l| ) includes fast local oscillations through the dependence on the harmonics 
and slow variation (modulation) in amplitude taken into account by the £, r 
dependence of u^. Introducing now this expansion and the slow variables in 
( JO| ) we may proceed to collect and solve different order 7 and I. We obtain: 
M o = u l = u 2 = u l = °7 v = l/k 2 , u i = a (C,t) 7^ and u 2 , = -4fc 3 /a. At 
order 7 = 3, I = 1 we obtain NLS for r) 

iA T + £-A^ + A\A\ 2 = 0. (5.2) 

The nature of solutions of NLS depends drastically of the sign of the product 
between the coeficient of A^ and that of A|A| 2 . In this case the product of 
a/k 3 by 4fc 3 /a is positive, and according to a well known stability criterion (see 
for exemple ||) the plane wave solution ( |2.4| ) is unstable. 



6 Conclusion and final coments on short waves. 

We have shown that ripples result from balance between linear dispersion and 
nonlinearity as in the case of long waves. However the physical model under 
consideration represents an ideal fluid because viscosity has been neglected, and, 
since dissipative phenomena take place at small scales, viscosity must affect 
asymptotic dynamics of ripples. 

The same approach can be used to study a flow of multiple layers. Such a 
flow occurs in many practical applications (e.g multilayer coating) in which the 
short wave behavior is of primordial importance. 

Note that, while the original Green-Naghdi system is Galilean invariant, 
this invariance is lost in the Camassa-Holm equation derived by Hamiltonian 
methods. In our case the perturbative theory conserves Galilean invariance and 
consequently equation ( |2.l| ) inherits this property from modified Green-Naghdi. 

In nonlinear dispersive systems, short wave dynamics does not occur as natu- 
rally as long wave dynamics. It is interesting to point out that some intermediate 
long wave shallow water models behave well, paradoxically, in the short wave 
limit. In these intermediate models a second asymptotic limit is always possi- 
ble, in general of long wave type, leading then to the ubiquitous KdV equation. 
In some cases however, there exists a second asymptotic limit of short wave 
type, leading to new nonlinear evolution equations. This is the case for BBMP, 
the integrable Camassa-Holm equation, and one of the Boussinesq system. The 
short wave limit of BBMP reads u xt = u — 3it 2 , where u has the same 

meaning as in fl2.ip . Its solitary wave solution comes from the balance between 
dispersion and nonlinearity as in KdV. Explicit solutions, blow-up and nonlin- 
ear instabilities are studied in ]l6| . The Camassa-Holm equation has as short 
wave limit, an integrable equation which belongs to the Harry-Dym hierarchy 
(k is a constant): u xt — ku — ^(ux) 2 — uu xx . Finally the short wave limit of the 
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Boussinesq system results as : u xU — u t — uu x . The linear limit for all these 
cases is the equation u x t — au (a constant) or else 



Ut — a udx = 0. (6-1) 

J — oo 

This linear nonlocal equation for u corresponds to the linear local wave equation 
u t — au x = that appears in the long wave case. Note that the short wave limits 
of BBMP, Camassa-Holm and Boussinesq system are not evolution equations 
stricto sensu, rather being integro-differential equations. Thus the nonlocality 
of nonlinea r ev olution equations for short waves appears to reflect the basic non 
locality of 
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